Quantum chaos in the mesoscopic device for the Josephson flux qubit 
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We study the quantum spectra and eigenfunctions of the three-junction SQUID device designed 
for the Josephson flux qubit at high energies. We analyze the spectral statistics on the parameter 
region where the system has a mixed classical phase space where regular and chaotic orbits can 
be found at the same classical energy. We perform a numerical calculation of eigenvalues and 
eigenstates for different values of the ratio of the Josephson and charging energies, Ej/Ec, which is 
directly related to an effective h parameter. We find that the nearest neighbour distributions P(s) 
of the energy level spacings are well fitted by the Berry-Robnik theory employing as free parameters 
the pure classical measures of the chaotic and regular regions of phase space in the different energy 
regions in the semiclassical case. The phase space representation of the wave functions is obtained 
via the Husimi distributions and the localization of the states on classical structures is analyzed. We 
discuss for which values of Ej /Ec it can be possible to perform experiments that could be sensitive 
to the structure of a mixed classical phase space. 

PACS numbers: 74.50.+r, 05.45.Mt, 85.25.Cp, 03.67.Lx 



I. INTRODUCTION 

In recent years, several types of superconducting qubits 
have been experimentally proposed r^&2£ These systems 
consist on mesoscopic Josephson devices and they are 
promising candidates to be used for the design of qubits 
for quantum computational 3 -! 4 -^! 7 - Indeed, a large ef- 
fort is devoted to succeed in the coherent manipulation of 
their quantum states in a controlable way. The progress 
made in this case allows to have nowadays Josephson 
circuits with small dissipation and large decoherence 
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Very recently, it has been proposed that, due 
to these developments, it could also be possible to use 
mesoscopic Josephson devices for the study of the quan- 
tum signatures of classically chaotic systems^! 9 - In Ref. H 
the quantum dynamics of the Device for the Josephson 
Flux Qubit (DJFQ) has been studied. In particular, it 
has been discussed how the fidelity (or Loschmidt echo)^ 
of the DJFQ could be studied experimentally for energies 
corresponding to the hard chaos regime in the classical 
limit. Here, we extend the work of Ref. by analyz- 
ing the possibility of studying, in the DJFQ, the mixed 
chaos regime (i.e., the energy range where there is a co- 
existence of chaotic and regular orbits in the classical 
limit). To this end, standard tools of analysis of "quan- 
tum chaos", like spectral statisticsil^^^i&^&i 9 . 
and phase space distributions ) 20 ' 21 ' 22 ' 23 will be used. 

It is by now well established that from the analy- 
sis of the spectral properties of quantum systems in 
the semiclassical regime it is possible to obtain infor- 
mation about the underlying dynamics of the classical 
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The probability distri- 



bution P(s) of the spacings s between successive energy 
levels - the nearest neighbor spacing distribution P(s)- 
unveils information on the associated classical dynam- 
ics. For integrable systems the levels are uncorrelated, 
and P(s) obeys a Poisson distribution. For completely 
chaotic classical motion, P(s) follows the prediction of 



the Random Matrix Theory (RMT)^i and when time 
reversal symmetry is preserved P(s) is closely approxi- 
mated by the Wigner distribution for the Gaussian Or- 
thogonal Ensemble (GOE), P(s) ~ s exp (— s 2 )^ 

Generic quantum systems do not conform to the above 
special cases, the classical phase space typically presents 
mixed dynamics, with coexistence of regular orbits and 
chaotic motioni 13 i 14 ' 15 i 16 ' 17 i 18 ' 19 In this generic case Berry 
and Robnik 13 proposed an analytical expression for the 
corresponding P(s), based on the knowledge of pure clas- 
sical quantities related to the Liouville measure of the 
chaotic and regular classical regions. The idea behind 
their calculations is that each regular or irregular phase 
space region gives rise to its own sequence of energy lev- 
els. For each region the level density results proportional 
to the Liouville measure of the classical region and the as- 
sociated level spacing distribution follows the Poisson or 
the Wigner form for regular and chaotic regions respec- 
tively. In the semiclassical limit these sequences of energy 
levels can be supposed independent and the complete dis- 
tribution P(s) is obtained by their random superposition. 
Several works have studied numerically the level statis- 
tics in systems with mixed dynamics ! 14 ' 15 ' 16 ' 17 ' 18 ' 19 Sys- 
tems with two degrees of freedom have been analyzed 
by several groups, mostly quartic oscillator a 14 ' 15 ' 16 and 
billiards ; 17 ' 18 ' 19 and in some works the Berry-Robnik pro- 
posal has been tested in detaih 16 ' 17 ' 18 

In contrast to the level statistics, the wave functions 
of quantum chaotic systems have remained relatively 
less explored. In particular the analysis of wave func- 
tions in phase space representations, such as the Wigner 
function^ or the Husimi distribution ) 21 ' 22 ' 23 allows a di- 
rect comparison between the classical and the quantum 
dynamics. Of particular interest are the zeros of the 
Husimi distribution which seem to be organized along 
regular lines or fill space regions for regular or chaotic 
classical dynamics respectivelyi 2 ^ 

Besides the importance of visualizing the dynamical 
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properties of quantum systems in phase space, techniques 
for measuring these functions, referred as "quantum to- 
mography " 25 ' 26 are subjects of active research in many 
experimental systems, like ion traps, optical lattices, en- 
tangled photons ; 27 ' 28 and also superconducting qubits. 29 

Josephson junctions have been used for the study of 
classical chaos since the early 1980s . 30 ' 31 A single under- 
damped junction with a periodic current drive can be- 
come chaotic in a wide range of parameter values^ - Sev- 
eral experiments have indeed studied this problem and 
measured chaotic properties in current-voltage curves 
and in voltage noise in Josephson junctions. 31 Moreover, 
networks with several junctions have been proposed for 
the study of spatio-temporal chaos. 32 All this cases cor- 
respond to classical chaos in dissipative systems with a 
time-periodic drive. Much less studied has been the case 
of classical hamiltonian chaos in Josephson junctions, 33 
mainly due to the fact that dissipation through a shunt 
resistance and/or coupling to the external measuring cir- 
cuitry is typically important. For the same reason, i.e., 
the difficulty in fabricating Josephson circuits with negli- 
gible coupling to the environment, few examples of quan- 
tum chaos in Josephson systems are found in the litera- 
ture. One of them is the work of Graham et al.^ who 
considered dynamical localization and level repulsion in 
a single Josephson junction with a time periodic drive. 
More recently T. D. Clark, M. J. Everitt and coworkers 3 ^ 
explored chaos and the quantum behaviour of SQUID 
rings coupled to electromagnetic field modes. The recent 
development of Josephson devices for quantum compu- 
tation, which need large coherence times, lead to signif- 
icant advances in the fabrication of circuits with small 
coupling to the external circuit and negligible dissipa- 
tion. This opened the possibility of using this type of 
mesoscopic devices for the study of quantum chaos. For 
example, Montangero et al& have proposed recently a 
Josephson nanocircuit as a realization of the quantum 
kicked rotator. The difficulty in realizing experimentally 
their system resides in that it needs to move mechani- 
cally one superconducting node in a high-frequency pe- 
riodic motion. A different proposal has been put for- 
ward in Ref|9j, where it has been shown that the Device 
for the Josephson Flux Qubit (DJFQ)^^! 7 - which con- 
sists on a three-junction SQUID, is classically chaotic 
at high energies. It could be therefore possible to use 
this system for the experimental study of quantum sig- 
natures of classical chaos. One possibility is the analy- 
sis of the fidelity or Loschmidt echo^ - in the quantum 
dynamics. 4 An experimental setup for the measurement 
of the Loschmidt echo in the DJFQ has been proposed 
in Ref|9|. In the above mentioned work, the system is 
prepared initially with a wave packelj 3 ^ localized in coor- 
dinate (phase) and momentum (charge) with an energy 
corresponding to the regime of hard chaos in the classical 
limit. The quantum evolution of the wave packet is evalu- 
ated in the unperturbed and the perturbed hamiltonians, 
and the overlap of the two evolved wave functions defines 
the Loschmidt echo or fidelity 10 , which can be measured 



experimentally. 9 Different behavior could be observed if 
the wave packet is initially localized in a chaotic or in a 
regular region of the phase space. Therefore, an interest- 
ing case to analyze is when the wave packet is prepared 
initially with an energy within the regime where there is a 
mixed phase space in the classical limit. In this case, one 
would expect that the behavior of the Loschmidt echo 
could depend on the location of the average coordinate 
and momentum of the initial wave packet. For exam- 
ple, in Refl37l it has been found a strong dependence of 
the fidelity with the initial state for mixed dynamics in 
the phase space in the case of Bose-Einstein condensates. 
However, in order to be sensitive to the structure of phase 
space in the case of mixed dynamics, it is necessary to 
have a small effective H. The aim of the present work is to 
analyze the quantum spectra and wave functions of the 
DJFQ in order to obtain for which values of the effective 
ft the quantum physics of this system can show signa- 
tures of the structure of the phase space in the case of 
mixed dynamics. To this end, we will use standard tools 
of quantum chaos theory by calculating numerically the 
level statistics of the DJFQ for different effective ft and 
the Husimi distribution. 

Concerning the spectral analysis, the quantum signa- 
tures of chaos have been discussed through the P(s) dis- 
tribution in Ref. [H for a SQUID with three junctions in 
the hard chaos regime. However, the case with only on- 
site capacitances was considered there (the capacitance 
of the junctions was neglected). Nevertheless, the de- 
vice for the Josephson flux qubit fabricated by the Delft 
goup£ has small on-site capacitances, about two orders 
of magnitude smaller than the intrinsic capacitances of 
the junctions^ This fact turns the model hamiltonian 
for the DJFQ to be different from the one studied in 
RefsIUll One of the goals of this paper is to analyze 
the spectral properties of the DJFQ considering realistic 
values of the different capacitances to analyze the device 
for the Josephson flux qubit (DJFQ) in the case of mixed 
dynamics. In addition we analyze the structure of the 
Husimi functions for the DJFQ, an issue that has been 
so far unexplored. The paper is organized as follows. In 
Sec|H]we introduce the quantum model for the device for 
the Josephson flux qubit. Before presenting the quantum 
spectral analysis, we will study in Sec. Illll thc dynamics 
of its classical analog. The presence of chaos will be char- 
acterized through the analysis of a measure of the chaotic 
volume, that will be defined and obtained as a function of 
the energy. We devote the rest of Sec. EI] to the analysis 
of the spectral properties. The NNS distribution will be 
obtained for different energies corresponding to different 
classical energy regions and dynamics and for different 
values of the effective H. In Sec. IIVI we compute the 
Husimi distribution for the DJFQ in order to character- 
ize the localization of the quantum states on typical phase 
space structures related to the different classical regimes. 
Finally in Sec. [V] we summarize our results and discuss 
possible experimental characterizations of the quantum 
manifestations of chaos in this system. 
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II. MODEL FOR THE DEVICE FOR THE 
JOSEPHSON FLUX QUBIT 

The DJFQ consists of three Josephson junctions in a 
superconducting ring 2 - that encloses a magnetic flux $ = 
/<f>o, with $o = h/2e, see FigHJ 




FIG. 1: Circuit for the Device for the Josepshon Flux Qubit 
as described in the text. Josepshon junctions 1 and 2 have 
Josepshon energy Ej and capacitance C, and junction 3 has 
Josepshon energy and capacitance a times smaller. The ar- 
rows indicate the sign convention for defining the gauge in- 
variant phase differences. The circuit encloses a magnetic flux 



The junctions have gauge invariant phase differences 
defined as tpi, cp2 and <£>3, respectivily, with the sign con- 
vention corresponding to the directions indicated by the 
arrows in Fig[TJ Typically the circuit inductance can 
be neglected and the phase difference of the third junc- 
tion is: if3 = —tpi + f 2 ~ 27r/. Therefore the system 
can be described with two dynamical variables: Lp\,<{>2- 
The circuits that are used for the Josephson flux qubit 
have two of the junctions with the same coupling energy 
Ej.i = Ej^ = Ej, and capacitance, C\ = C2 = C, while 
the third junction has smaller coupling -Ej.3 = aEj and 
capacitance C3 = aC, with 0.5 < a < 1. The above 
considerations lead to the Hamiltonian 2 -^ 



n 



1 



(1) 



where the two-dimensional coordinate is tp = (^1,^2)- 
The potential energy is given by the Josephson energy of 
the circuit and, in units of Ej, is: 

V(ip) — 2+a~cos(pi — cos(p 2 — acos(2Trf+(pi—ip2) ■ (2) 

The kinetic energy term is given by the electrostatic en- 
ergy of the circuit, where the two-dimensional momen- 
tum is 



P=(Pi,P 2 



M 



dCp 
~dt' 



and M is an effective mass tensor determined by the ca- 
pacitances of the circuit, 



M 



V27T 



with 



m 



1 + a + 7 —a 
—a 1 + a + 7 



We included in M the on-site capacitance C g = jC. 
(Typically 7 ~ 10~ 2 — 1CP 3 1). In the presence of 
gate charges Q 9t i induced in the islands, the momentum 
is P -> P + ^Q g 2 The system modelled with Eqs. flTJ- 
@ is analogous to a particle with anysotropic mass M 
in a two-dimensional periodic potential V(0)M- 

In typical junctions, the Josephson energy scale, Ej, 
is much larger than the electrostatic energy of electrons, 
Ec = e 2 /2C, and the system is in a classical regime. On 
the other hand, mesoscopic junctions (with small area) 
can have Ej ~ Ec, and quantum fluctuations become 
important^ 2 - In this case, the quantum momentum oper- 
ator is defined as 

P -» P = -ihv v = -ih(—, —). 

oipi aip2 



After replacing the above defined operator P in the 
Hamiltonian of Eq.(fT]), the eigenvalue Schrodinger equa- 
tion becomes 



V(0) 



9(0) = E9{0) , (3) 



where we normalized energy by Ej and momentum by 
h/yJ&Ec/Ej. We see in Eq.([3]) that the parameter 
7] = \J 8Ec I Ej plays the role of an effective h. It is 
well-known that the ratio Eq/Ej controls the effect of 
quantum fluctuations in single Josephson junction a 43 i 44 
and in arrays of several Josephson junctions. 45 ' 46 For 
Ej ^> Ec, {rj 1), the junctions can be described 
with a classical dynamics; while for Ej ~ Ec, (v ~ 1) 
the effect of quantum fluctuations becomes important. 43 
Experiments where the Josephson junctions are fabri- 
cated for different values of Ec /Ej have been performed 
both for single junctions^ 4 - and for junction arrays.— In 
the last case quantum phase transitions as a function of 
Ec/Ej have been studied. 45 ' 46 Therefore, the parameter 
rj = y/8Ec/Ej is a natural choice for quantifying the 
effective h in this system. 

For quantum computation implementations 2 -' 6 -'^ the 
DJFQ is operated at magnetic fields near the half-flux 
quantum (/ = fo + df, with /o = 1/2). For values of 
a > 1/2, the potential Eq.([2]) has two well defined min- 
ima. At the optimal operation point / = 1/2, the two 
lowest (degenerated) energy states are symmetric and 
antisymmetric superpositions of two states correspond- 
ing to macroscopic persistent currents of opposite sign. 
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The offset value Sf determines the level splitting between 
these two states. These eigenstates are energetically sep- 
arated from the others (for small Sf) and therefore the 
DJFQ has been used as a qubit^&i (i.e. a two-level trun- 
cation of the Hilbert space is performed). In addition 
the barrier for quantum tunneling between the states de- 
pends strongly on value of a and its height goes up as a 
is increased. The possibility to manipulate the potential 
landscape by changing a is a crucial point for experi- 
mental implementation of qubits. Typical experiments 
in DJFQ have values of a in the range 0.6 — 

As we will discuss here, the higher energy states of the 
DJFQ show quantum manifestations of classical chaos. 
In what follows we focus our study of the DJFQ consid- 
ering the realistic case of: (i) small on-site capacitances, 
taking 7 = 0.02, (ii) a magnetic field corresponding to 
the optimal operation point of the DJFQ, / = 1/2, and 
(iii) the values of a = 0.7 and 0.8 in coincidence with the 
experimental values employed in Ref. l6j|7|. 

III. SPECTRAL STATISTICS 



0.2- (a) 




Before entering into the analysis of the quantum spec- 
tra we will focus on the classical dynamics of the DJFQ. 
As we already anticipated in the Introduction, generic 
systems present mixed classical dynamics and the DJFQ 
is not the exception. Therefore for a given energy E our 
aim is to estimate the chaotic volume v c -t(E), defined as 
the probability of having a chaotic orbit (i.e. Lyapunov 
exponent A > 0) at energy E. As we will show below, 
this parameter will be relevant in the statistical analysis 
of the quantum spectrum. 

The classical dynamical evolution was obtained solving 
the Hamilton equations derived from Eq. (fT]) : 

m.^ = -V^(£), (4) 

where we have normalized energy by Ej and time by t c — 
a/ h 2 C/4:e 2 Ej = hjr\Ej (the Josephson plasma frequency 
is ojp — t^ 1 ). The numerical integration was performed 
with a second-order leap-frog algorithm with time step 
At = 0.02i c . 

For different values of the parameter a and magnetic 
flux / we compute the maximum Lyapunov exponent A 
for each classical orbit at different energies E. We esti- 
mate the chaotic volume V(±(E) using 10 3 initial condi- 
tions chosen randomly with uniform probability within 
the available phase space for each given energy. Also 
the average Lyapunov exponent, X(E), of the chaotic or- 
bits is obtained. These results are shown in Fig. [2] for 
a = 0.8 and / = 1/2. We observe that both v c h(E) and 
X(E) increase smoothly with energy, as it is usual in sev- 
eral similar systems with two degrees of freedomi 47 ' 48 i 49 
Above the minimum energy of the potential, E m i n , we 
find: (i) regular orbits for E m i n < E < E c h (w c h = 0), 
(ii) soft chaos (i.e., coexistence of regular and chaotic or- 
bits, < D c h < 1) for E c h < E < Ehc with the average 



FIG. 2: (a) Average maximum Lyapunov exponent A and (b) 
chaotic volume v c h versus energy E for a = 0.8 and / = 1/2. 



Lyapunov exponent A > above E c h and (iii) hard chaos 
(all orbits are chaotic, v c h — 1) for E > E^. The bound- 
aries of these different dynamic regimes as a function of 
a, in the range [0.5, 1.0_L and /, in the range [0, 0.5], have 
been obtained in Ref. 0- Here we will focus on the case 
with / = 1/2 and we will study some different cases of 
a. 

In order to look for signatures of quantum chaos, we 
follow a standard statistical analysis of the energy spec- 
trum. First we calculate the exact spectrum {E n } by 
diagonalization of the quantum Hamiltonian. The eigen- 
value equation Eq. ([3]) is solved by discretizing the phases 
with Aip = 27r/1000, and the resulting hamiltonian ma- 
trices of size 10 6 x 10 6 are diagonalized using standard 
algorithms for sparse matrices. We have verified that in- 
creasing the discretization by a factor of 2 does not affect 
the results of the spectrum within the needed accuracy 
for the ranges of energies studied here. As we mentioned 
we set 7 = 0.02 and / = 1/2, and we obtain eigenvalue 
spectra for different values of the parameters r\ and a 
defined in the previous section. 

The level spectrum is used to obtain the smoothed 
counting function N av (E) which gives the cumulative 
number of states below an energy E. In order to anal- 
ize the structure of the level fluctuations properties one 
unfolds the spectrum by applying the well kwown trans- 
formation x n — N av (E n )*£ From the unfolded spectrum 
one can calculate the nearest-neighbor spacing (NNS) 
distribution P(s), where = Xi+i — Xi is the NNS. 



We have taken into account the symmetries of the 
Hamiltonian Eq.([T]). For / = 1/2 the Hamiltonian has 
reflection symmetry against the axis ip2 = fi and against 
the axis ip2 — —<p\- The eigenstates can be chosen with 
a given parity with respect to these two symmetry lines. 
Therefore, we compute the NNS distribution employing 
eigenstates of a given parity. This kind of decomposition 
is a standard procedure followed in the analysis of spec- 
tral properties of quantum systems whenever the Hamil- 
tonian of the system possesses a discrete symmetry^ We 
consider the even-even parity states and the NNS distri- 
bution is computed for different energy regions inside the 
classical interval (E c h, Ehc), corresponding to soft chaos, 
and for energies E > E\ lc ( and E < 2A), corresponding 
to hard chaos. 

The Berry- Robnik theory seems to be suitable to an- 
alyze, in the semiclassical regime, sequence of levels of 
quantum systems whose classical analogous presents co- 
existcncy of regular and chaotic dynamics (i.e., soft chaos 
regime) . If p\ and P2 are the relative measures of the reg- 
ular and chaotic parts of the classical phase space then, 
the Berry-Robnik distribution^ 3 - reads: 

P BR (s) = p\ exp {-pis) erfc Qv^p 2 s^ + 

(2P1P2 + ^ W P2 S J ex P (-Pis - ^pls 2 ^j , (5) 

where pi + p% = 1. It is easy to verify that P BR (s) 
interpolates between the Poisson and Wigner GOE dis- 
tributions as — > p\ — > 1, but does not exhibit level 
repulsion for p\ ^ 0. 

In Fig[3] we show the cumulative level spacing distri- 
bution W(s) = j P(s)ds obtained numerically following 
the prescription described before. We have done this in 
order to describe in some detail the behavior for small 
values of s, (in the following we denote the cumulative 
distributions by the same name that the corresponding 
NNS distribution) . In all the cases we have fitted the nu- 
merically obtained W(s) employing Eq.© for the NNS 
distribution, and we have extracted the fitted quantum 
parameter p\ = pb r - 

The particular results presented in Fig. [3] correspond 
to a window of ~ 100 eigenvalues around E c h < E — 
1.6 < Ehc, within the soft chaos regime, Fig. [3] (a),(b); 
and Ehc < E — 2, within the hard chaos regime, Fig. [3] 
(c),(d). We take the realistic experimental value for the 
parameter a = 0.8 and consider different values of the 
quantum parameter rj : the case with rj — 0.01 is shown 
in Fig. [3] (a),(c); and the case with 77 = 0.05 is shown 
in Fig. [3] (b),(d). We should remark that the classical 
dynamics is independent of the parameter 77, which has 
a pure quantum origin and plays the role of an effec- 
tive Planck's constant in the Schrodinger equation, as 
we mentioned before. In addition in Fig[3J we show for 
comparison the W(s) corresponding to the Poisson and 
Wigner GOE distributions. 

We first discuss the case with 77 = 0.05, that is al- 
ready smaller than in the cases studied in Refill where 
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FIG. 3: Cumulative distribution W(s) for a — 0.8 and 
/ = 0.5. See the text for details. The continuous line is the 
fitted Berry-Robnik distribution. We show for comparison 
the Poisson (dotted line) and Wigner (dashed line) cumula- 
tive distributions. Top panels correspond to E = 1.6 with 
(a) rj — 0.01 and (b) rj = 0.05. Bottom panels correspond 
to E = 2.0 with (c) rj = 0.01 and (d) rj = 0.05. The fitted 
Berry-Robnik parameters are (a) pb r = 0.44, (b) pb r = 0.93, 
(c) p br = 0.99 and (d) p br = 0.96. 



ij = 0.07 — 0.17 was considered. The numerical results 
for rj = 0.05, in the case of hard chaos [E — 2.0, shown 
in Fig. [3J (d)], are in good agreement with the Wigner 
distribution, and we obtain pt r — 0.96. In a case corre- 
sponding to mixed classical dynamics [E = 1.6, shown in 
Fig. E](b)], we find that the distribution departs slightly 
from the pure Wigner form. However, we have obtained 
Pbr — 0.93 ^> v c h ~ 0.4, meaning that the level distri- 
bution in this case does not seem to be very sensitive 
to the mixed phase space expected in the classical limit. 
The reason is that for increasing 77 the mean energy level 
spacing increases (proportional to rj 2 for large energies), 
and therefore the width of the energy region evaluated for 
the statistics with a given number of levels (~ 100 in this 
case) also increases in the same way. Since v c h (E) varies 
rapidly with E within the soft chaos region, relating its 
value with the fitted pb r , which is obtained evaluating 
the statistics over a wide energy region, becomes mean- 
ingless for large 77. Indeed, deep in the quantum regime 
the Berry-Robnik fitted parameters are not expected to 
be related to the classical measure of the chaotic (regular) 
part of the phase space. 13 ' 16 ! 17 

We now discuss a smaller value of the effective H, corre- 
sponding to rj — 0.01. In Fig. [3](a), for E — 1.6 (mixed 
classical dynamics), we find now that the W(s) clearly 
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FIG. 4: Fitted Berry-Robnik parameter pf, r as a function of 
the dimensionless energy E/Ej for / = 0.5 and rj = 0.01. 
(a) a = 0.7 and (b) a = 0.8. The horizontal error bars in 
energy are defined by the interval of eigenenergies used in 
the statistics, and it is a decreasing function of the density 
of states. The vertical error bars correspond to error in the 
parameter pf, r as obtained from the numerical fits. The dotted 
line shows the chaotic fraction of the classical phase space v c h 
obtained from the classical dynamics. 



departs from the pure Wigner form, and that it can be fit- 
ted with the Berry- Robnik distribution with pb r = 0.44. 
This value is very close to the classical chaotic volume 
for this case, v c h ~ 0.4. In the case for E = 2 (hard 
chaos), shown in FigJ3] (c), we have obtained pb r = 0.99, 
in agreement with v c h — 1 and also in good agreement 
with the Wigner distribution, as expected i 12 ' 14 

In general we find that in a nearly semiclassical regime, 
r\ = 0.01, the numerical results for the Berry-Robnik pa- 
rameter pi, r show a good agreement with the classical 
measure p\ , that by definition is equivalent to the chaotic 
volume v c h- This is analyzed in Fig. 3] where we plot the 
quantum parameter pb r obtained for different sections of 
the spectra with ~ 100 eigenvalues around a given en- 
ergy E. We show results for two cases of the parameter 
a. The chaotic fraction of the classical phase space v c h is 
also plotted. The results for pb r and v c h are very close to 
each other. When changing the parameter a the location 
in energy of the onsets of the soft chaos and hard chaos 
regimes shifts. We also see that the curves of pb r vs. E 



shift in the same way, giving further support to the cor- 
respondence between p br and v c h ' 14]16]17]18 These results 
corroborate the validity of the Berry-Robnik theory in 
the semiclassical energy region that corresponds to small 
effective Planck's constant, as it is the case for r\ = 0.01. 

Besides the cases reported above, we have also ana- 
lyzed a few other values of a in the range 0.5 — 0.9 and 
/ = 0.4, 0.5, obtaining similar results for the spectral 
statistics. In general, we observe that in order to obtain 
a spectral statistics with a Berry-Robnik paremeter that 
agrees with the classical measure of the chaotic volume 
values of 77 < 0.05 are needed. 



IV. PHASE SPACE AND HUSIMI 
DISTRIBUTIONS FOR THE DJFQ 

In this section we pursue our study of the signatures 
of quantum chaos presenting an analysis of the quantum 
phase-space distributions in the case of mixed classical 
dynamics. Taking into account the analysis of the previ- 
ous section we focus on the case 77 = 0.01. 

Quantum phase space distributions are of increasing 
interest in studies of quantum chaos because they allow 
a direct comparison between classical and quantum dy- 
namics. The Husimi distribution associated to a quan- 
tum wave function \^>) (see definition below, Eq.( [6])) 
it is based on the coherent-state representation and is 
well suited to represent wave functions in phase space 
because it is always real and possitivei 21 i 22 i 23 ' 24 Due to 
these properties it is usually referred as a quasi probabil- 
ity distribution. 

In order to compute the Husimi function for the DJFQ 
we must take into account the fact that the classical 
phase space is four dimensional. The Husimi distribu- 
tion function for a state \^f) is 



CPo,V>bl*) 



(G) 



where |Poj<A)) corresponds to minimum- uncertainty 2tt- 
periodical wave packets^ given by 



\Po,<Po) 



C x exp[iK ■ (0 — <p )] x 

cos((/? ,i - ipi) + cos(v? ,2 



(7) 



exp 



<P2) 



2a 2 



where Kq = (ki,k2) with fcijfe integers and Po = r/Ko. 
The width of the wave packet is a = \frjjs, with s the 
squeezing parameter, and we choose the value s = 3.23, 
which is the same value used in Ref(9j for the initial co- 
herent wave packets. 

The potential has two minima for / = 1/2 which are at 
(<p*,—(p*) and (—<p*,<p*), with cosip* = l/2a. To better 
analyze the Husimi function, it is convenient to make the 
following change of variables: 



fx = 



V2 • 



Pr 



P, 



Pl-P2 

P1 + P2 

V2 • 



(8) 
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In this way the two minima lie along the direction of ip x . 
The normalization by y/2 is chosen such that new vari- 
ables satisfy [ip x ,P x ] — irj, [(p y ,P v ) — if] in the quantum 
regime. 

The classical Poincare surface of section is calculated 
in the plane (cp x ,P x ), taking ip y = and P y > 0. We 

want to compare the Husimi distribution p^j (K, (p) corre- 
sponding to the eigenstate \^ v ) with eigenvalue E v with 
the classical Poincare section at an energy E « E v . To 
this end, we construct an analog of the surface of sec- 
tion by obtaining a two-dimensional section of {K , ip) 
(which is a four-dimensional density in phase space) in 
the following way^ 



<f>X(P x ,<p x )=p*(P x ,pZ;cp x ,0) 



(9) 



0, Pjf is ob- 



where, given the values P x ,(p x and <p y 
tained such that the classical energy is equal to E and 
the possitive root, Py > 0, is chosen. 

We obtain numerically the eigenstates \^f v ) of Eq. ([3]), 
after using a discretization of Aip = 27r/500. Then, us- 
ing Eqs. ([6])- ([9]), we compute the sections of the Husimi 
distributions, (P x ,(p x ). In order to characterize the 
localization of the quantum states on the classical phase 
space structures, we choose a few examples of for 
eigenstates that lie in energy regions corresponding to 
regular classical dynamics E < E c h and soft chaos re- 
gion, E ch < E < E hc , respectively. 



show here three cases corresponding to eigenstates with 
energies E v = 1.5219, E v = 1.5208 and E v = 1.5193 
(panels (a) , (b) and (c) respectively). The localization 
of these states on the stability islands and fixed points is 
clearly observed. 
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FIG. 6: Classical Poincare surface of section for E = 1.6. 
Section of Husimi phase space distribution, $^ (P x , ip x ) for 
(b) E v = 1.601, (c) E v = 1.6008, (d) E v = 1.5993. 



1.0 




FIG. 5: (a) Classical Poincare surface of section for E — 
1.52. Sections are symmetric with respect to tp x — > — ip x and 
P x — » — P x ; only the region of ip x > and P x > is shown. 
Section of Husimi phase space distribution, ^ (P x , (p x ) for 
(b) E v = 1.5219, (c) E v = 1.5208, (d) E v = 1.5193. 

In Fig. [5] (a) we plot for E = 1.52 < E c h the classical 
Poincare section in which the stability islands associated 
to the regular dynamics are observed. We have computed 
the Husimi phase space distributions $^ (P x , ip x ) for sev- 
eral eigenstates (~ 20) near the energy E — 1.52. We 



In FigE^a) and FigJTJa) we plot for classical ener- 
gies E — 1.6 and E — 1.7 respectively, the clas- 
sical Poincare sections together with a selection of 
some of the calculated Husimi phase space distribu- 
tions ^{P x , (p x ) for eigenstates with energies E v = 
1.601,1.6008,1.5993 (FigE](b), (c) and (d) respectively) 
and E v = 1.6992,1.7004,1.6994 (FigH (b), (c) and (d) 
respectively) . 

In these cases the soft chaos behavior is evident by 
the structure of the Poincare sections, in which regular 
islands are sorrounded by chaotic regions. The local- 
ization of the states on classical structures like already 
distroyed chains of islands is observed in the figures. In 
addition, the Husimi distribution of Figj7](d) corresponds 
to a state localized on the chaotic region of Fig[7]Ja). 

The above analysis of the Husimi distributions shows 
that for rj — 0.01, it is possible to use localized wave 
packets as initial conditions for the experimental mea- 
surement of the Loschmidt echo;&2& since they can sense 
the structure of the phase space with mixed classical dy- 
namics in this case^i 



V. CONCLUSIONS 

In this paper we have characterized the quantum sig- 
natures of chaos in the three-junction SQUID device. For 
realistic parameter values the classical dynamics exhibits 
different regimes that go from mixed dynamics to fully 
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FIG. 7: Classical Poincare surface of section for E = 1.7. 
Section of Husimi phase space distribution, ^ (P x , tp x ) for 
(b) E v = 1.6992, (c) E v = 1.7004, (d) E v = 1.6994. 



developed chaotic motion. As a consequence the spectral 
statistics, characterized by the distribution of the nearest 
neigbour energy spacing (NNS) P(s) in the high energy 
region, is expected to unveil signatures of the mentioned 
behavior. The analysis has been performed for different 
energy regions inside the classical intervals corresponding 
both to the soft chaos (i.e., mixed phase space) and hard 
chaos regimes, and we considered the even-even parity 
states to compute the NNS distribution. Our numerical 
results show that, for rj < 0.05 (and for r\ = 0.01 in partic- 
ular), in a nearly semiclassical regime, P(s) is well fitted 
by Berry-Robnik like formulae, where the pure classical 
measures of the chaotic and regular regions have been 
used as the only free parameters. 

We also found that the individual eigenstates can also 
be intimately linked to the phase space structures that 
characterizes the different classical regimes for rj < 0.05. 
In order to analyze how quantum states are supported or 
localized on different classical structures that are present 
in the different regimes in this case, we have investigated 
the Husimi phase space distributions for different eigen- 
states with energies E v in the classical interval. We would 
like to mention that there are few studies of Husimi dis- 
tributions for Hamiltonian systems with two degrees of 
freedom,22, as it is the case of the DFJQ studied in the 
present work. 

One important advantage of Josephson junction de- 
vices is that they can be fabricated with well-controlled 

parameters. The effective h, is h c g — rj = i and 



since Ej oc A and Ec oc 1/A, with A the area of the 
junctions, we have that H e s oc 1/A. Thus, the fabrica- 
tion of different DJFQ with junctions with varying area 
could allow to study cases with H e s spanning from the 



semiclassical to the quantum regime. This is indeed im- 
portant since different regimes can be accessed experi- 
mentally depending on the magnitude of r\. The qubit 
regime of two-level dynamics of the DJFQ is observed 
experimentally in devices with r\ m 0.4j&1 In Ref[§ it has 
been found that signatures of chaos in the Loschmidt echo 
can be observed at high energies E ~ 3Ej in devices with 
an effective h of the order of r\ ~ 0.1. Here we have shown 
that the observation of the quantum effects in the case 
with mixed chaotic and regular orbits (for an intermedi- 
ate energy range) needs the study of devices in a more 
semiclassical regime with r\ w 0.01. This could motivate 
experimental measurements looking for the dependence 
of the Loschmidt echo^ with initial conditions, due to the 
phase space structure of the mixed classical dynamics, if 
the experiments are performed in devices with rj « 0.01. 
Considering the values^ of Ej ~ 250GHz ~ 2K and the 
operation temperature of 20 mK reported by the Delft 
group£, typical level spacings of O.OlEj ~ 20 mK can 
be experimentally resolved in the device of Ref@. This 
energy resolution is enough for the case of the Loschmidt 
echo in devices with 77 ~ 0.1, analyzed originally in Ref0. 
However, the semiclassical regime explored in this work 
(77 = 0.01) requires a resolution in the level spacings of 
the order of 5 x 1Q~ 4 Ej. Thus, for experiments in the 
cryogenic range (20 mK) devices with larger values of Ej 
should be employed. On one hand, a smaller 77 ~ 0.01 
already requires junctions with larger area A, and there- 
fore larger Ej. On the other hand, Josephson junctions 
fabricated with high T c superconductors^ can have a 
large Ej. Therefore, devices designed with high T c super- 
conductors can be good candidates for the experimental 
challenge of studying the mixed phase space in the semi- 
classical regime of the DJFQ. 

Another possible type of experiment is to start the 
system in the ground state and apply a constant pulse 
in some external parameter (for instance, the magnetic 
field). After the pulse is applied, the probability of re- 
maining in the ground state could be related to the en- 
ergy level statistics^ Also, an interesting experiment 
could be to perform studies of the low frequency noise, as 
it has been done in mesoscopic chaotic cavities i 53 ' 54 For 
example, one could drive the DJFQ into the hard chaos 
regime with a voltage source such that Ey = \CV 2 > 
Ehc (and V < 2A/e) and then measure the noise in the 
current. How the current noise is related to the spectral 
statistics in this case is a very interesting problem, which 
could be the subject of future studies. 
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